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Using the quark-meson coupling (QMC) model, we study nuclear matter from the point 
of view of quark degrees of freedom. Performing a re-definition of the scalar field in matter, 
we transform QMC to a QHD-type model with a non-linear scalar potential. The potentials 
obtained from QMC are compared with those of the relativistic mean-field models. 



Since at present rigorous studies of quantum chromodynamics (QCD) are limited 
to matter system with high temperature and zero baryon density, it is important to 
build models which help to bridge the discrepancy between nuclear phenomenology 
!>■ \ and QCD. We have proposed a relativistic quark model for nuclear matter and finite 

nuclei, that consists of non-overlapping nucleon bags bound by the self-consistent 
■ exchange of isoscalar, scalar (cr) and vector {cu) mesons in mean-field approximation 

I (MFA) - this model is called the quark-meson coupling (QMC) model. On the 

' other hand, recent theoretical studies show that various properties of finite nuclei 

. can be very well described by the relativistic mean-field (RMF) models, i.e., Quantum 

hadrodynamics (QHD).cP In this letter, we consider relationship between QMC and 
QHD and study how the internal structure of the nucleon sheds its effect on effective 
^ . nuclear models. 

^ I In our previous works, the MIT bag model has been used to describe the quark 

C ■ structure of the nucleon. Since the confined quarks interact with the scalar field, a, 

in matter, the effective nucleon mass, M'^, in QMC is given by a function of a 
^ I through the quark model of the nucleon. (Although the quarks also interact with 

\^ • the u! nieson, it has no effect on the nucleon structure except for a shift in the nucleon 

energy. B*) The (relativistic) constituent quark model (CQM) is an alternative model 
for the nucleon. Recently, Shen and TokiB^ have proposed a new version of QMC - 
the quark mean-field (QMF) model, where CQM is used to describe the nucleon. 

In the present study, as well as the bag model (BM), we want to use the rela- 
tivistic CQM with confining potentials, V{r), of a square well (SW) and a harmonic 
oscillator (HO) to see the dependence of the matter properties on the quark model. 
It is assumed that the light (u or d) quark mass, rriq, is 300 MeV in CQM, while 
ruq = MeV in BM. Furthermore, we introduce a Lorentz-vector type confining 
potential, which is proportional to 70, as well as the scalar one: 

V{r) = {l+P^o)U{r), (1) 

where the potential, U{r), is given by SW or HO and (3(0 < /3 < 1) is a parameter to 
control the strength of the Lorentz-vector type potential. We assume that the shape 
of the Lorentz-vector type confining potential is the same as that of the scalar type 
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one. 

In SW, the solution for a quark field, V'gi can be calculated d la Bogolioubov.a) 
The potential is given by U{r) = for r < R and M for r > R, where R is the 
radius of the spherical well and M is the height of the potential outside the well. 
After finishing all calculations, the limit M ^ oo is taken. eP This system may be 
described by Lagrangian density 



Csw = Mn ■ d - m,)ij,e{R - r) - + /?7 • a)Vg5(r - R), (2) 



where a'^ is the unit vector in time direction: = (1,0). This Lagrangian provides 
a boundary condition: i'y ■ nipg = (1 + (3^ ■ a)ipq at r = R, where n'^ is the unit normal 
outward from the potential surface. This condition gives 



, , (l-P)(E-mn) , , „ 

with jn the spherical Bessel function, x the eigenvalue of the confined quark and E 
the quark energy. (Note that in fact the quark current fiows in only the azimuthal 
direction. Il*) To take into account the corrections of the spurious center of mass (cm.) 
motion and gluon fiuctuations, we add the-faxniliar form, —z/R, with a parameter 
z for those corrections to the total energy. □^'B' In SW, the nucleon mass, Mn, (at 
rest) is then given by M^r = (3a — z)/R, where = + and = Rmg. 

In case of HO, we set /? = 1 in Eq.dl|) because the quark wave function can be 
obtained analytically. 0^ When U{r) = ^cr^ (c the oscillator strength), a condition 
to determine the quark energy is given by 



E + mq{E - mg) = 3^/c. (4) 



le non-relativistic harmonic oscil- 
) Thus, the nucleoli mass is given 



The cm. energy can be evaluated exactly, as in 
lator, and it is just one third of the total energy. 

by Mn = 2E — Eg, where Eg describes gluon fiuctuation corrections. 0^ 

For the MIT bag model there are many good reviews. 0^ In BM, we take = 

MeV and (3 = 0. (Even in BM it is possible to include the Lorentz-vector type 

potential using Eqs.(|2|) and (^. However, if we use a large (3 in BM, it is hard to 

get good values of the nuclear matter properties.) 

Now we consider an iso-symmetric nuclear matter with Fermi momentum kp, 

which is given by ps = 2kp/3Tr'^ {pB the rmclear matter density). Then, the total 

energy per nucleon, Etot, can be written ascP 

Etot = 3 / " dk^M*^\a)+P + + ^ps, (5) 

where is calculated by the quark model. The o" and u meson masses, rrio- and m^, 
are taken to be 550 MeV and 783 MeV, respectively. The a; field is determined by 
baryon number conservation: uj = QuPB/fnt idi^ is the w-nucleon coupling constant), 
while the scalar mean-field is given by a self-consistency condition: (dEtot/da) = O.IIP 
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In SW, we set the radius of the potential to be i? = 0.8 fm and detrmine z so 
as to fit the free nucleon mass, Mjv (= 939 MeV). The parameter f3 is chosen to be 
and 0.5 to examine the effect of the Lorentz-vector type confining potential. We 
find that z = 4.396 and 5.164 for /3 = and 0.5, respectively. In HO, there are two 
adjustable parameters, c and Eg. We determine those parameters so as to fit the free 
nucleon mass and the root-mean-square (charge) radius of the free proton: r% = 0.6 
fm^.cP {rN is calculated by the quark wave function.) We find that c = 1.591 fm~^ 
and Eg = 344.7 MeV for the free nucleon. In nuclear matter, we keep c and Eg 
constant and the quark energy E varies, depending on the scalar field. In BM, the 
bag constant, B, and the parameter, z, are fixed to reproduce the free nucleon mass. 
As in SW, we choose the bae radius of the free nucleon to be 0.8 fm. We find i?^/^ 
= 170.3 MeV and z = 3.273.11) 

Now we are in a position to de- 
Table I. Coupling constants, and K. The termine the coupling constants: the a- 
effective nucleon mass, M^, is calculated nucleon coupling constant, g^, and gl 

at On. The nuclear incompressibility, K, is n a j. cj. J.^ i t_- j- 

" . ^, , , . ' ' are hxed to nt the nuclear bmdmg en- 

quoted in MeV. The SW model with f3 = ^ ■ ^ 

0(0.5) is denoted by SW0(5). ^rgy (-15.7 MeV) at the saturation den- 

sity (po = 0.15 fm ) for nuclear matter. 
The coupling constants and some calcu- 
lated properties for matter are listed in 
Table I. The present quark models can 
provide good values of the nuclear in- 
compressibility, K. 

In SW and BM with massless quarks, the quark scalar density in the nucleonll* 
vanishes in the limit /5 — > 1, which means that the a meson does not couple to the 
nucleon. EP This fact implies that as f3 is larger the u-nucleon coupling is weaker in 
matter. Thus, we can conclude that qualitatively a large mixture of the Lorentz- 
vector type confining potential leads to a weak scalar mean-field and hence a large 
effective nucleon mass in nuclear matter. Since in MFA a small effective nucleon 
mass (and hence a strong scalar field) is favorable to fit various properties of finite 
nuclei, I-* the confining potential including a strong Lorentz-vector type one may not 
be suitable for describing a nuclear system. 

The main difference between QMC and QHD at hadronic level El) lies in the 
dependence of the nucleon mass on the scalar field in matter. By performing a 
re-definition of the scalar field, the QMC Lagrangian densityB* can be cast into a 
form similar to a QHD-type mean-field model, in which the nucleoiL mass depends 
on the scalar field linearly, with self-interactions of the scalar field. t3) In QMC, the 
nucleon mass in matter is given by a function of a, Qj^jfj{a), through the quark 
model of the nucleon, while in QHD the mass depends on a scalar field linearly, 
^NQHD ~ — go4> {4' is the scalar field in a QHD-type model). Hence, to 
transform QMC into a QHD-type model, we can apply a re-definition of the scalar 
field, 

go(t){a) = Mn - M^QMci^), (6) 





gi 




MI,/Mn 
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swo 


84.4 


104 


0.725 


329 


SW5 


66.6 


65.2 


0.807 


287 


HO 


147 


64.5 


0.805 


309 


BM 


67.6 


66.1 


0.805 


278 



to QMC, where go is a constant chosen so as to normalize the scalar field (j) in the 
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limit (J ^ 0: (l>{a) = a + 0{(7'^). Thus, go is given by go = -{dM^Qjy.ifj/da)^=o- In 
QMC, we find 50 = for SW and BM, while go = ^g^ for HO. 

The contribution of the scalar field to the total energy, E^ci, is now rewritten in 
terms of the new field d) 



Escl = 11 df [(Va)2 + mla^] = J 



df 



/i(0)2(V0)2 + C/,(0) 



where [/, describes the self-interactions of the scalar field 



[/,((/,) = lmla{(t>)^ and h{(f)) 



dUs{(t>) 



(7) 



(8) 



Note that in uniformly distributed nuclear matter the derivative term in Egci does 
not contribute^The effect of this term on the properties of finite nuclei has been 
studied in Ref.ll9).) Now QMC can be re-formulated in terms of the new scalar field, 
(f), and it is of the same form as QHD with the non- linear scalar potential, Us{(j)), 
and the coupling, /i((/>), to the gradient of the scalar field. (Note that since this 
re-definition of the scalar field does not concern the vector interaction, the energy of 
the uj field (see Eq.(D) is not modified.) 

The Zimanyi-Moszkowski (ZM) modelM-' is a good example. By re-definig the 
scalar field, ZM can be exactly transformed to a QHD-type niodel with a non-linear 
potential. Since the effective nucleon mass in ZM is given byEJ) 



Mn 



l + {g^a/MN)'' 



(9) 



the model involves higher order couplings between the a and the nucleon. Introduc- 
ing a new scalar field cp by go4>{(^) = — Afjv ZA/(^)' easily find 



90 = 9a, (t>{cr) = ^ , and (t(0) = - 
1 + qo 1 



with g' = gcr/M^. Thus, the non-linear potential is given by Eq.(^) 

2 



(10) 



(11) 



In general, the in-medium nucleon mass may be given by a complicated function 
of the scalar field. However, in QMC the mass can be parametrized by a simple 
expression up to 0{g1):a 



M^/Mn ^l-ay + by\ 



(12) 



with y(= gf^a/Mj\[) a dimensionless scale and two (dimensionless) parameters, a and 
b. This parametrization is accurate up to ~ 4/9o-EJ' 

Once the parameters, a and 6, are fixed, we can easily re-define the scalar field 
using Eq.(P). We find 



go = aga, 4>{^) = c — da"^ and cr((/)) 



1 - yi - Adcj) 
2d ' 



(13) 
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with d = bgu/aM]^. This satisfies the condition: a — > in the hmit 
non-hnear potential is thus calculated 



0. The 



U, 



2 



d 



^1.2, / 'n^l 



Mn) 



+ 0{gl) 



(14) 



where r = b/a. 



Table II. Parameters a, 6, k and A. 



The standard form of the non-linear 
scalar potential is usually given by 





a 


b 


K(fm-^) 


A 


swo 


1.01 


0.215 


19.2 


78.8 


SW5 


1.02 


0.497 


39.0 


327 


HO 


0.687 


0.245 


42.3 


384 


BM 


0.998 


0.435 


35.1 


264 
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more sophisticated version of QHD, 
inspired by modern methods of ef- 
fective field theory, many other terms of 
meson-meson and meson-nucleon couplings are considered. However, we here focus 
on only the parameters k and A to make our discussion simple. It is well known that 
the non- linear scalar potential, Eq.(|l5|), is practically indispensable to reproduce the 
bulk properties of finite nuclei and nuclear matter in RMF. 

We can estimate the parameters, k and A, in QMC by comparing Eqs.(14) and 
(15). In Table || the parameters, a, b, k and A, are presented. The QMC model 
leads to the non-linear potential with k ~ 20 — 40 (fm~^l and A ~ 80 — 400. Since 
QMC predicts that both a and b are always positive,^'^ we can expect that the 
quark substructure of the in-medium nucleon provides a non-linear potential with 
positive K and positive A in the QHD-type mean-field model. 

Using those parameters, we can draw the non-linear scalar potential generated 
by QMC, which is illustrated in Fig. 1. The non-linear pDtentials, which have been 
used in various relativistic mean-field models, lliP'll3'll3)'EJ)'ll3) are also presented in 
Fig. 2. From Fig. 1 we can see that the quark models lead to the similar non- linear 
potentials, in spite of the big difference in the confinement mechanism. On the 
other hand, the non-linear potentials in RMF show quite different behaviors. The 
potentials in NLBI13) and PMSEj) are relatively close to those produced by QMC. 

In summary, we have calculated the properties of nuclear matter using QMC 
with various quark models for the nucleon. Then, we have performed a re-definition 
of the scalar field in matter and transformed QMC to a QHD-type model with a non- 
linear scalar potential. The QMC model gives k ~ 20 — 40 (fm~^) and A ~ 80 — 400 
for the non-linear scalar potential. The shapes of the potentials generated by the 
quark models are very similar to one another, although the confinement mechanism 
is quite different in each model. On the contrary, the parameters, n and A, phe- 
nomenologically determined in RMF take various values and the potentials in RMF 
are quite different from one another for large \(p\. In general, the phenomenological 
potential may consist of the part, which is caused by the quark substructure of the 
nucleon, and inherent self-couplings of the scalar field in matter. It is very intriguing 
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Fig. 1. Non- linear scalar potentials generated 
by QMC. The solid curve shows Us ~ 
^^(j)^. The dashed curve (with solid cir- 
cles) is for SW with /3 = 0(0.5), while the 
result of BM is shown by the dotted curve. 
The dot-dashed curve is for HO. 
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Fig. 2. Non-linear scalar potentials in RMF. 
The solid curve shows Us ~ -^0^- The 
dashed curve, with open (solid) circles is for 
G2 (NLB),y' while the dashed one with- 
out any marks is for PM3.li3 rlihe dashed 
curve with crosses is for ZM-tif The po- 
tentials in TMlEJ and NLltJ are respec- 
tively shown by the dotted and dot-dashed 
curves. 



if the potential due to the internal structure of the nucleon could be inferred by 
analyzing experimental data in the future. 
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